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Thick  GaN  crystals  can  now  be  grown  (by  the  non-equilibrium  CVD  reaction 
of  GaCl  and.  NH£)  thar  are  free  of  the  precipitates,  voids  and  ancrocracks 
that  plagued  us  last  year-  However*  they  are  still  dominated  electrically 
by  a native  shallow  donor  effect.  A high  pressure-high  temperature  system 
for  annealing  crystals  and  for  growing  crystals  by  liquid  phase  epitaxy 
under  high  pressure  Nh  instead  of  Nttt  has  been  assembled  and  calibrated, 
and  the  first  experiments  begun  to  efiminate  the  native  donor,  so  as  to^  (over) 
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Introduction 


GaN  Is  a wide,  direct  gap  III-V  semiconductor,  with  a bandgap  of  3.4eV 
In  the  ultraviolet,  and  which  may  prove  to  be  a useful  semiconductor  for 
microwave  and  for  electroluminescence  devices.  Its  possible  usefulness  for 
power  microwave  amplifiers  (high-field,  transit-time  controlled  devices)  is 
based  upon  a predicted  high  saturated  drift  velocity  and  upon  a predicted 

high  pair-production  threshold.  Taken  together,  these  predicted  values 

i 

yield  a flgure-of -merit  20  times  higher  than  that  for  Si  for  a fast  power 
amp! ifier. 

The  goal  of  this  contract  has  been  the  preparation  of  suitable  samples 
to  measure  these  two  vital  parameters,  and,  incidentally,  to  establish  the 
technological  base  for  the  use  of  GaN  In  semiconductor  devices.  However, 
three  problems  have  interfered  with  the  achievement  of  these  goals. 

First,  the  melting  point  of  GaN  is  in  the  neighborhood  of  2000°C  and 
the  equilibrium  N2  pressure  over  such  a melt  may  be  greater  than  40,000  atm. 
Even  at  1000°,  well  below  the  melting  point,  the  decomposition  pressure  is 
already  100  atm.  The  trick  to  growing  GaN  is  to  note  that  at  1000°,  TOO  atm 

of  N„  should  be  in  equilibrium  with  10  torr  of  NH,.  However,  the  homogeneous 

2 3 o 

dissociation  of  NH^  into  N2  and  H2  is  very  slow  at  temperatures  below  1100 

and  so  low  pressures  of  NH3  can  be  used  instead  of  high  pressures  of  N2  to 

grow  GaN  well  below  its  melting  point,  in  our  case,  by  its  chemical  vapor 

. 

P reaction  with  GaCl.  During  this  process,  the  GaN-  slowly  dissociates  into  Ga 

[and  N2-  Hence,  the  growth  is  basically  a non-equilibrium  process,  where  we 

force  the  formation  reaction  to  be  faster  than  the  decomposition  reaction. 
Second,  unless  the  crystals  are  intentionally  compensated,  their 


electrical  properties  are  dominated  by  the  presence  of  a native  shallow 
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donor  at  concentrations  between  10  and  10  per  cm  and  no  low  resistivity 
p-type  material  has  ever  been  reported. 

Third,  GaN  is  normally  grown  on  sapphire  substrates.  Thin  crystals, 
less  than  about  30um  thick,  often  exhibit  grain  boundaries,  voids,  occlusions, 
high  dislocation  densities  and  a considerable  degree  of  strain.  These  de- 
fects tend  to  diminish  as  the  thickness  increases.  In  addition,  the  resis- 


tivity tends  to  increase  with  increasing  thickness. 

In  previous  reports,  we  surmised-  that  the  native  donors  might  be  associ- 
ated in  some  way  with  these  defects,  and  so  we  described  the  rapid  growth  of 
thick  (up  to  5mn)  crystals.  However  these  were  characterized  by  very  heavy 
Ga  precipitation  and  by  the  appearance  of  cracks  and  voids.  Last  year  we 
described  a high  pressure-high  temperature  system  that  we  were  assembling  in 
order  to  anneal  GaN  crysals  and  to  grow  GaN  by  liquid  phase  epitaxy  using  N2 
instead  of  NH-j  at  temperatures  much  closer  to  the  true  melting  point.  The 
system  was  designed  to  go  to  1500°C  and  10,000  atm,  of  N2,  so  as  to  study 
and  control  the  native  donor,  assuming  its  density  to  be  relatable  to  equi- 
librating N2  pressure  and  assuming  1500°  to  be  high  enough  to  allow  equili- 
bration with  bulk  crystals.  At  that  time  we  noted  that  all  available  crys- 
tals were  degenerately  doped  n-type  by  the  native  donor,  so  that  the 
electrical  properties  due  to  intentionally  added  impurities  could  not  be 
determined  by  simple  Hall  and  resistivity  measurements  vs.  temperature. 

This  meant  that  the  basic  semiconductor  doping  parameters,  which  elements 
are  donors,  which  acceptors,  which  shallow,  which  deep,  were  not  known  for 
GaN.  However,  we  discovered  a technique,  low  temperature,  bound-exciton 
luminescence,  whichwasable  to  detect  at  least  two  acceptors  and  one 
shallow  donor  in  addition  to  the  native  shallow  donor,  unmasked  by  the 
pervasiveness  of  the  native  donor. 
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In  this  report  we  show  that  the  native  donor  is  not  associated  with 

crystal  imperfections.  We  can  now  grow  thick  crystals  on  sapphire  via  GaCl 

and  NHg  that  are  not  cracked  and  are  relatively  free  of  precipitates  and  voids. 

We  can  now  grow  thinner  crystals  (to  30um)  in  which  resistivity  and  defects 

no  longer  depend  on  thickness.  We  can  grow  extremely  planar  surfaces  If  we 

use  other  substrate  orientations.  However,  in  all  of  these,  the  onrl- 

present  native  donor  still  dominates  the  electrical  properties. 

% 

The  high  pressure-high  temperature  system  is  now  set  up  and  calibrated 

to  1200°  and  5,000  atm.  of  N2-  We  have  succeeded  in  growing  GaN  by  liquid 

phase  epitaxy  grom  Ga  and  N2  and  have  started  to  anneal  previously  grown 
specimens  at  these  high  N2  pressures.  It  is  still  too  early  to  tell  whether 
we  can  affect  the  native  donor  by  these  treatments. 

We  have  doped  CVO  grown  GaN  crystals  by  ion  implantation  or  by  doping 
during  growth  with  a wide  selection  of  elements  predicted  to  be  shallow 
donors  and  shallow  acceptors,  or  to  be  expected  as  common  impurities.  This 
was  done  partly  to  establish  a semiconductor  technology  for  GaN  - what  are 
the  shallow  donors  and  acceptors  - and  partly  in  an  attempt  to  compensate 
or  otherwise  affect  the  native  donor.  The  low  temperature  photo luminescence 
experiments  extended  to  these  materials  have  now  provided  us  with  a wealth 
of  information  about  these  doping  elements,  not  only  from  bound-exci ton 
luminescence,  but  also  from  donor-acceptor  and  free-to-bound  transitions. 

We  have  now  determined  for  the  first  time  which  elements  are  shallow  donors 
and  which  shallow  acceptors  In  GaN, what  their  ionization  energies  are,  and 
also  which  elements  give  rise  to  deep  levels. 

Ouring  the  photo! uminescence  experiments,  we  discovered  a broad  emission 
band  which  dominated  the  spectra  at  the  highest  excitation  intensities. 

This  band  exhibited  stimulated  emission  gain.  By  analyzing  this  band  (spectra, 
excitation  dependence  of  spontaneous  and  stimulated  emission,  dependence  on 
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excitation  volume,  temperature  dependence)  we  were  able  to  show  that  the 
emission  was  due  to  formation  of  electron-hole  liquid  metallic  droplets.  We 
were  able  to  work  out  the  phase  diagram  of  this  electron-hole  liquid  and  to 
show  that  this  phase  Is  stable  to  over  80°K,  the  highest  stability  of  such 
a phase  in  any  material.  By  preparing  good  optical  cavities  in  several 
crystals  we  were  able  to  make  ultraviolet  lasers  via  this  recombination 
mechanism  using  optical  pumping.  Finally,  from  threshold  and  gain-loss 
measurements  on  such  lasers,  we  were  able  to  show  that  the  losses  deter- 
mined indicate  that  these  GaN  crystals  are  of  good  semiconductor  quality, 
despite  the  presence  of  the  native  donor. 

The  net  result  of  our  progress  this  year  is  the  realization  that  it 
is  only  the  native  donor  that  stands  in  the  way  of  our  making  definitive 
drift  velocity  and  pair  production  measurements  and  of  utilizing  GaN  as  a 
"normal"  semiconductor.  But  now  perhaps,  we  have  the  tools  to  alter  the 
native  donor  density  (high  pressure  ^ at  high  temperatures)  and  the  means 
to  measure  it  (low  temperature  photoluminescence). 

Standard  Crystal  Growth  by  Chemical  Vapor  Deposition 

Our  previous  reports  have  described  the  slowly  improving  quality  of 
GaN  grown  epitaxially  on  sapphire  by  the  CVD  reaction  of  GaCl  and  NH^  at 
about  1050°C.  Our  early  results  Indicated  that  the  most  rapid  growth 
occurred  on  R-plane  ( Toi 2 ) substrates,  although  the  surface  growth  morpho- 
logy was  not  planar.  We  learned  how  to  clean  the  sapphire  (Hj  at  1200°)  to 

avoid  an  extra  slow-growth  stage  required  during  nucleation.  We  learned  how 
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to  minimize  contamination  in  the  10  -10  cm  range  due  to  Si,  a shallow 
donor,  by  using  an  alumina  liner  in  the  furnace  to  eliminate  contact  with 
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the  fused  silica  envelope  and  by  using  higher  purity  starting  materials.  The 

electrical  properties  were  now  dominated  by  an  always-present  shallow  donor, 
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2-5  x 10  cm  , which  we  showed  to  be  a native  defect  by  eliminating  any  associ- 
ation between  it  and  any  expected  contaminant  in  that  concentration  range. 

Since  the’  bulk  resistivity  of  the  crystals,  increased  with  thickness,  we  surmised 
that  this  donor  might  be  due  to  growth  related  imperfections  which  decreased  with 
thickness,  and  we  began  a program  to  grow  very  thick  crystals.  By  optimizing  the 
temperature  profile,  the  flow  rates,  substrate  placement  and  flow  geometry,  we 
were  able  to  grow  GaN  single  crystals  on  2x2cm  R-plane  sapphire  up  to  10mm 
thick.  Last  year  we  set  up  a second  CVD  apparatus  to  grow  such  thick  crystals 
routinely.  Last  year  we  also  reported  that  such  rapidly  grown  thick  crystals 

were  characterized  by  three  types  of  defects:  (1)  a pervasive  precipitation  of 
o 

very  small  (<  100  A)  Ga  particles,  (2)  large  voids  frequently  containing  free  Ga 

and  usually  oriented  on  crystallographic  planes,  and  (3)  microcracks  on  crystal 

cleavage  planes,  starting  from  the  uppermost  surface,  but  seldom  passing  through 

the  entire  thickness  of  the  epitaxial  GaN  layer.  These  defects  lead  to  catas- 
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trophic  breakdown  even  in  very  modest  electric  fields  (10  -10  V/cm),  and 
aborted  our  attempted  measurements  of  the  high  field  drift  velocity. 

This  year,  we  have  therefore  devoted  a major  effort  in  an  attempt  to 
understand  the  origin  of  these  defects  in  order  to  eliminate  them.  Our  tools 
have  been  optical  microscopy,  scanning  electron  microscopy,  transmission  elec- 
tron microscopy  and  back-reflection  Laue  x-ray  diffraction. 

First  we  consider  the  cracks.  The  coefficient  of  thermal  expansion  of 
sapphire  is  larger  than  that  of  GaN.  Thus  on  cooling  the  newly  grown  crys- 
tal from  1050°C  to  room  temperature,  the  sapphire  shrinks  more  than  the  grown 
crystal  and  therefore,  the  sapphire  should  be  in  tension,  the  GaN  in 
compression.  Hence  it  is  the  sapphire  that  should  crack,  not  the  GaN. 
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In  fact,  the  sapphire  does  crack.  Whenever  the  GaN  layer  thickness  is 
greater  than  about  70um,  the  sapphire  substrates  always  exhibit  cracking 
and  crazing  extending  completely  through  their  thickness.  It  is  observed 
that  the  GaN  does  not  crack  at  all  unless  its  thickness  is  at  least  200um 
and  then  the  microcracks  formed  radiate  from  the  top  surface,  rather  than 
from  the  interface  with  the  sapphire  and  these  usually  do  not  propagate 
all  the  way  through  to  the  interface.  A very  simple  and  straight-forward 
stress  analysis  shows  why  the  GaN  cracks  even  though  it  should  be  under 
compressive  stress.  The  differential  strain  arising  at  the  Interface  re- 
sults in  a bending  moment  with  the  outer  surface  of  the  GaN  being  bowed 
into  a convex  form.  Although  the  stress  in  the  GaN  is  compressive  at  the 
interface,  the  bowing  causes  it  to  slowly  change  to  a tensile  stress  at  the 
outer  surface  and  this  stress,  maximum  at  the  outer  surface,  increases  for 
thicker  cyrstals  where  the  bending  (bowing)  moment  increases.  This  then 
qualitatively  explains  why  the  thick  layers  (>  200pm)  exhibit  microcracks 
and  why  these  originate  at  the  outer  surface,  rather  than  at  the  interface. 
Quantitatively,  the  calculated  minimum  tensile  force  which  causes  cracks  to 
appear  at  the  onset  200pm  thickness  due  to  differential  thermal  contraction 
of  the  GaN-sapphire  sandwich  from  1050°  to  room  temperature  is  about  two 
orders  of  magnitude  less  than  the  predicted  ultimate  fracture  strength  of 
an  ideally  perfect  GaN  crystal.  Thus  the  cracks  should  originate  at 
local! ized  crystallographic  defects.  In  fact,  this  is  exactly  what  is 
observed.  The  cracks  all  originate  at  the  upper  GaN  surface  and  always 
pass  through  one  or  more  of  the  large  voids  (discussed  below)  lying  on  or 
just  below  the  surface.  Thus,  to  grow  crystals  thicker  than  200pm,  growth 
should  be  interrupted  prior  to  attaining  this  thickness.  The  sapphire  sub- 
strate should  then  be  removed  (chipped  and  ground  off)  and  growth  resumed 
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without  the  presence  of  the  GaN-sapphire  interface.  In  fact,  this  procedure 
does  work  and  very  thick  crystals  can  be  grown  without  the  associated 
microcracks.  However,  as  noted  below^  the  original  motivation  for  growing 
very  thick  crystals  (the  increase  in  resistivity  with  thickness)  is  no 
longer  operative.  For  relatively  pure  crystals,  primarily  low  Si  contami- 
nation, the  resistivity  no  longer  increases  with  thickness,  as  long  as  the 
thickness  is  greater  than  30-50um. 

The  second  problem  associated  with  the  rapid  growth  of  thick  GaN  crys- 
tals on  R-plane  sapphire  is  the  frequent  occurrence  of  large  voids,  often 
containing  some  free  Ga.  We  have  now  shown  that  this  feature  is  associated 
both  with  the  substrate  orientation  and  with  the  fast  rate  of  growth.  The 
best  lattice  match  between  GaN  and  R-plane  (1012)  sapphire  is  the  (1120)  prism 
plane  of  GaN  and  the  GaN  does  exhibit  that  orientation,  with  its  c-axis  in 
the  plane  of  growth  oriented  parallel  with  the  projection  of  the  sapphire  c- 
axis  in  that  plane.  This  GaN  prism  plane  is.  a very  slow  growth  plane,  and 
unfortunately,  the  fast  growth  planes  are  not  perpendicular  to  it.  The 
latter  lead  to  the  development  of  other  (slow-growing  or  "equilibrium")  facets 
which  are  predominantely  R-planes  (1012)  or  R-like  planes  (lOln)  with 
smaller  angles  between  the  plane  and  the  c-axis-  This  results  in  a non-planar 
surface  morphology.  The  surface  consists  of  a series  of  ridges  bounded  by 
these  R-like  planes.  During  the  initial  growth  stage,  many  nuclei  form, 
each  bounded  by  such  planes.  As  the  nuclei  grow  and  merge,  these  form  a 
very  large  number  of  small  ridges  and  valleys  and  the  many  re-entrant 
boundaries  in  the  valleys  allow  rapid  growth  to  occur.  As  growth  progresses, 
the  ridges  tend  to  become  fewer  In  number  presumably  due  to  slightly  different 
growth  rates  in  adjacent  valleys.  As  the  ridges  merge,  the  R-like  facets 
become  larger  and  the  number  of  re-entrant  edges  decrease,  the  latter  re- 
sulting in  a decreased  over-all  growth  rate.  In  addition,  prism  plane  (slow 

growth)  facets  begin  developing  from  the  edges  of  the  crystals,  again 
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decreasing  the  number  of  re-entrant  edges  and  decreasing  the  growth  rate. 
Thus,  for  thick  crystals,  the  growth  mechanism  changes.  The  slow  growing 
R-like  planes  become  the  dominant  growth  planes.  These  planes  grow  with 
different  growth  velocities.  Thus,  for  a ridge  which  is  the  intersection  of 
two  different  planes  of  this  type,  one  side  of  the  ridge  grows  faster  than 
the  other  side.  This  leads  to  the  development  of  "overhanging  cliffs". 
Continued  growth  extends  these  overhangs  into  cave-like  structures,  which 
sometimes  get  covered  over,  forming  cavities.  These  are  the  large  gross 
voids  we  otre.'ve.  Once  the  cavi  ty  is  formed,  it  is  no  longer  in  contact  with 
NH^  and  the  GaN  slowly  decomposes  into  Ga  and  leaving  free  Ga  in  the 
voids.  This  problem  then  is  associated  with  the  appearance  of  well  developed 
large  area  facets  as  the  crystals  become  very  thick.  For  thick  crystals  the 
problem  cannot  be  el  iminated,  but  i t can  be  minimized  by  slowing  the  growth 
rate  down  and  supplying  excess  NHg  as  the  crystals  get  thicker.  This  slows 
the  growth  of  the  overhangs  and  allows  the  valley  growth  to  keep  up  with 
ridge  growth  and  avoid  the  formation  of  caves.  At  over-all  growth  rates  of 
60um/hr,  this  becomes  a problem  when  the  crystals  are  only  30um  thick.  At 
30um/hr  void-free  crystals  can  be  grown  to  60um  thicknesses,  and  for  very 
slow  growth  rates,  1-lOum/hr,  crystals  at  least  200um  thick  can  be  grown 
without  such  voids.  Again,  note  that  we  no  longer  desire  to  grow  crystals 
thicker  than  about  60um. 

The  last  major  defect  of  our  standard  CVD  grown  crystals  on  R-plane 
sapphire  is  the  appearance  of  a mist  of  very  small  Ga  precipitates  ( <1 00  R) 
scattered  throughout  the  volume  of  the  crystal.  These  scatter  light  and 
render  the  thicker  samples  completely  opaque  to  transmitted  light.  These 
defects  are  also  associated  with  growth  rate.  In  our  CVD  system  the 
growth  rate  is  directly  related  to  the  GaCl  flow  rate  (HC1  over  Ga  at  950°) 
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net  tothe  NH^  How  rate.  For  rapid  growth  of  very  thick  crystals  we 
increased  the  Ga  transport  rate,  not  the  NHg  flow.  Hence,  it  is  not  sur- 
prising that  under  such  conditions  free  Ga  forms  on  the  growing  surface  and 
gets  incorporated  into  the  growing  crystal  before  it  can  react  with  NH^.  We 
should  be  able  to  avoid  this  disproportionation  by  increasing  the  NH2  flow 
rate,  but  in  the  apparatus  we  now  have  we  are  already  working  close  to  its 
maximum  through-put  rate  for  NH^.  Hence  we  should  decrease  the  Ga 
flow  rate  and  thus  the  growth  rate.  This  does  indeed  work.  Crystals 
grown  at  growth  rates  between  1 and  20um/hr  are  free  of  such  precipitates. 
Note  that  such  lowered  growth  rates  are  also  needed  to  reduce  the  voids 
discussed  above. 

Thus  we  are  now  able  to  grow  very  thick  GaN  layers,  greater  than 
200um,  on  R-plane  sapphire.  A two-step  process  (removing  the  sapphire) 
eliminates  the  microcracks  and  slow  growth  eliminates  the  precipitation 
and  greatly  reduces  the  number  of  voids.  These  are  now  good  single  crys- 
tals, without  grain  boundaries  or  other  gross  defects.  The  only  structural 
defects  that  can  be  observed  are  microscopic  basal  plane  stacking  faults. 

Since  Si  was  removed  as  the  major  electrically-active  impurity  two 
years  ago,  the  resistivity  of  the  GaN  no  longer  increases  with  thickness 
once  the  samples  are  10-30\im  thick  and  it  is  only  the  native  donor  which 
controls  the  carrier  density.  Thus  we  no  longer  have  the  impetus  to  grow 
very  thick  crystals,  even  though  we  are  now  able  to  do  so.  Our  second  CVD 
apparatus  for  the  routine  growth  of  GaN  on  R-plane  sapphire  now  produces 
such  uniform  crystals  consistently  at  a growth  rate  of  20y/hr. 

Finally,  we  note  that  the  removal  of  all  the  defects  discussed  above 

had  no  effect  on  the  presence  of  the  native  donor.  The  native  donor  density 
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in  the  routinely  grown  crystals  is  10  -10  cm  and  the  mobility  is 

2 

consistently  about  100  cm  /V-sec. 


CVO  Growth  Experiments:  Basal  Plane  Growth,  Carbon  Doping 


Two  years  have  now  passed  since  we  made  the  decision  to  concentrate 
on  CVD  growth  on  R-plane  sapphire.  Growth  on  that  sapphire  orientation  yielded 
the  fastest  growing  crystals.  Since  then,  we  have  learned  how  to  best  clean 
the  substrate  surface,  how  to  eliminate  Si  contamination,  how  to  avoid 
cracks,  voids  and  precipitates  and  the  optimum  conditions  of  temperature, 
flow  rates  and  sample  placement  geometry,  to  produce  good  crystals.  With 
this  information,  and  knowing  that  many  of  our  R-plane  growth  problems  were 
associated  with  the  resultant  non-planar  surface  growth  morphology,  we  de- 
cided to  re- investigate  CVO  growth  on  other  sapphire  orientations.  Avoiding 
the  conditions  that  normally  result  in  precipitation  and  microcracks,  we 
found  that  of  several  orientations  tried,  basal  plane  sapphire  gave  excel- 
lent results.  Basal  plane  GaN  grows  on  basal  plane  sapphire.  Nucleation 
forms  islands  with  basal  plane  surfaces.  These  eventually  merge  when 
these  islands  are  at  different  heights,  producing  a planar,  but  stepped, 
surface  morphology.  Good  growth  proceeds  and  the  morphology  does  not 
change.  Growth  rates  of  30-60ym/hr  produce  good  crystals  as  defect  free 
as  the  best  R-plane  grown  crystals  and  the  morphology  is  planar.  The 
planar  surfaces, free  of  voids  and  overhangs, allow  us  to  avoid  grinding  and 
etching  steps  for  preparing  samples  for  electrical  measurements.  The  native 
donor  density  (10^-10^cm"^)  and  mobility  is  identical  with  those  in  R- 
plane  grown  samples.  Thus  we  will  now  be  growing  basal  plane  GaN  as  well 
as  GaN  on  R-plane  sapphire  in  our  standard  CVO  growth  apparatus.  So  far,  the 
only  disadvantage  of  this  is  that  basal  plane  sapphire  substrates  cost 
three  times  as  much  as  R-plane  wafers. 

Last  year's  Progress  Report  described  our  initial  attempts  at  using  low 
temperature  near  gap  photoluminescence  to  characterize  impurities  in  as- 
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grown  GaN.  One  emission  line  was  tentatively  Identified  with  carbon  on  the 
basis  of  circumstantial  evidence:  carbon  should  be  an  expected  Impurity  and 
the  line  was  enhanced  whenever  graphite  was  used  in  the  growth  apparatus. 
(This  identification  is  confirmed  this  year  via  the  ion- implantation  re- 
sults given  below.)  The  photoluminescence  line  was  identified  as  a bound 
exciton  line,  but  its  energy  was  too  low  to  assign  it  as  either  a shallow 
donor  or  as  a shallow  acceptor.  Thus  it  had  to  be  relatively  deep.  From 
naive  chemical  arguments  (covalent  radius,  electronegativity)  we  suspected 
that, as  a simple,  isolated,  substitutional  defect,  carbon  should  appear  on 
a N site  rather  than  on  a Ga  site,  and  thus  act  as  an  acceptor.  We  were 
looking  for  simple,  easily  handled  acceptors,  to  compensate  the  native  donor 
always  present  in  our  crystals.  Carbon  had  the  advantage  that  it  was 
easily  compatible  with  our  CVD  growth  method.  Thus  we  decided  to  dope  our 
standard  crystals  during  growth  by  introducing  a metered  flow  of  CH^  into 
the  reactor  as  a source  of  carbon.  The  photoluminescence  of  these  samples 
showed  a strong  increase  in  the  line  already  tentatively  identified  as 
carbon,  showing  that  the  identification  was  correct.  However,  the  native 
donor  density  did  not  exhibit  any  gross  changes  as  determined  by  its  bound 
exciton  spectra,  nor  was  there  any  change  in  the  shallow  donor  concentra- 
tions, as  measured  by  the  free  electron  Hall  coefficient  at  room  temperature. 

However,  there  was  an  unexpected  increase  of  the  electron  mobility  from  about 
? 2 

100  cm'/V-sec  to  almost  300  cm  /V-sec.  From  the  bound  exciton  spectra,  car- 
bon cannot  be  a shallow  donor  or  a shallow  acceptor.  If  It  is  a deep  accep- 
tor it  should  either  compensate  the  native  donor,  leading  to  a lower  free 
electron  concentration  and  a lower  mobility,  or,  if  a chemical  mass  action 
equilibrium  is  operative  at  the  growth  temperature,  it  should  enhance  the 
shallow  donor  density,  but  decrease  the  mobility.  If  carbon  is  a deep 
donor  it  should  either  not  affect  the  native  donor  nor  the  room  temperature 
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mobility,  or  If  a chemical  mass  action  equilibrium  occurs,  it  should  de- 
crease the  native  donor  density  leading  to  a lower  electron  density  at 
room  temperature  but  a higher  mobility.  None  of  these  four  possibilities 
agrees  with  the  experiments.  It  is  as  if  carbon  were  neither  a donor  nor 
an  acceptor,  but  simply  gettered  the  crystal,  improving  the  mobility.  If 
so,  the  luminescence  line  could  not  be  due  to  a point  defect.  Rather  the 
carbon  would  have  to  be  part  of  an  extended  defect  - perhaps  if  it  is 
basically  an  acceptor,  paired  with  the  native  donor  to  yield  an  electrically 
neutral  isoelectronic  center.  Although  carbon  does  not  lead  to  lower 
carrier  densities,  its  role  as  an  enhancer  of  mobility  is  intriguing  and 
is  being  pursued  further. 

High  Pressure  Synthesis  and  Processing 

The  central  theme  of  our  program  now  is  aimed  at  eliminating  or  at 
least  reducing  the  native  shallow  donor  defect,  perhaps  a nitrogen  vacancy. 
Our  approach  is  to  grow  crystals  or  to  anneal  CVO  grown  crystals  under 
varying  pressures  of  N^  rather  than  NH^  and  at  temperatures  that  are  suffi- 
ciently high  to  allow  rapid  diffusion  and  equilibration.  Under  such  condi- 
tions, >1100°,  NH^  is  no  longer  stable. 

During  the  past  year  we  have  received  the  components  of  the  high 
pressure-high  temperature  system  that  was  ordered  last  year.  This  con- 
sists of  a high  pressure  gas  compressor  (14,000  atm),  two  internally  heated 
vessels  capable  of  sustaining  1400°C  at  10,000  atm  and  an  externally  heated 
vessel  (1100°  at  4,000  atm).  These  are  in  addition  to  a 3,000  atm  compressor 
we  already  had. 

The  system  was  installed  with  its  auxiliary  equipment,  tested  section 
by  section  and  then  together.  The  ixperiments  were  delayed,  originally  by 
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malfunctions  due  to  shipping  damage,  manufacturers'  errors  and  our  unfanrl- 
llarlty  with  this  equipment , and  later  by  problems  arising  out  of  contami- 
nation of  gas,  gas  leakage  in  the  vessels  and  in  the  compressor,  corrosion 
of  parts  due  to  contamination,  failure  of  the  compressor  due  to  gas  scouring 
of  the  compression  cylinders  and  other  miscellaneous  "bugs". 

In  the  first  experiments,  the  externally  heated  pressure  vessel  was 
used  to  anneal  CVO  grown  GaN  crystals  at  temperatures  between  800°  and  1030° 
and  M2  pressures  up  to  3,000  atm.  Care  was  exercised  to  remove  oxygen  and 
water  vapor  from  the  ambient.  These  experiments  had  two  purposes:  to 
deduce  an  accurage  pressure- temperature  phase  diagram  for  the  three-phase  GaN- 
Ga-N2  system  to  predict  conditions  for  high  pressure  synthesis  and  annealing, 
and  secondly,  to  begin  studying  the  direct  effects  of  high  pressure  ^ 
anneals. 

SEM  analysis  of  the  GaN  surfaces  before  and  after  annealing  provided 
no  evidence  of  whether  or  not  the  samples  were  dissociating,  since  similar 
features  appeared  on  the  surface  irrespective  of  whether  the  annealing  was 
done  well  above  or  well  below  the  estimated  equilibrium  pressure.  Low 
temperature  photoluminescence  experiments  showed  no  marked  change  in  the 
shallow  donor  density  due  to  annealing.  These  negative  results  in  these 
early  experiments  at  relatively  low  pressures  and  temperatures  are  not  un- 
expected. They  were  designed  to  give  us  experience  with  the  use  of  the 
equipment  and  with  the  surface  effects  on  GaN  crystals  while  the  equipment 
was  being  tested  and  calibrated  into  its  much  higher  pressure  and  temperature 
regimes. 

Recently,  experiments  were  begun  in  one  of  the  two  internally  heated 
vessels  at  higher  pressures,  taking  precautions  to  preclude  contamination 
by  O2  or  I^Q.  A fused  silica  crucible  containing  Ga  and  occasionally  a sapphire 
substrate  or  a CVD  grown  GaN  crystal  is  placed  in  the  furnace.  After 
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applying  pressure  (4,000  to  8,600  atm)  the  furnace  is  raised  to  the 
annealing  or  growth  temperature  (960-1180°)  held  there  for  3-8  hours,  then 
rapidly  cooled  and  depressurized.  The  Ga  melt  surface  was  invariably 
covered  with  a brittle  crust  containing  sharply  faceted  crystallites  (vlOym 
in  size)  and  a large  number  of  whiskers  as  seen  by  SEM  and  optical  microscopy. 
The  surfaces  of  the  sapphire  substrates  were  covered  by  a uniform  but  thin 
layer  (10-20um)  of  faceted  crystal.  The  areas  of  the  samples  were  too  small 
for  good  x-ray  identification  or  clear-cut  electrical  characterization,  but 
energy  dispersive  analysis  in  the  SEM  indicated  that  Ga  is  the  dominant  con- 
stituent (N  is  not  detectable)  with  Si  and  A1  as  trace  contaminants.  So 
far  the  electrical  measurements  indicate  high  resistivity  (which  is  what  we 
are  seeking),  but  on  these  thin  polycrystalline  layers  this  could  be  due  to 
the  grain  boundaries.  SEM  analysis  of  GaN  grown  on  CVD  GaN  substrates  did 
not  indicate  unequivocal  features  that  could  be  attributed  to  either  growth 
or  dissolution  of  the  substrate.  Therefore  separate  experiments  are  now 
being  done  to  characterize  the  surface  of  CVD  grown  substrates  etched  by 
dissolution  at  low  pressures  and  in  contact  with  Ga  to  catalyze  dissolution. 

Preparations  are  now  being  made  to  grow  GaN  from  a Ga  solution 
of  with  the  solution  in  a thermal  gradient. 


Ion  Implantation 
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GaN  is  Invariably  heavily  doped  n-type  with  10  -10  native  shallow 
donor  defects  per  cm3.  At  room  temperature  these  are  all  ionized  yielding 
a corresponding  concentration  of  degenerate  free  electrons.  Thus  the  Hall 
coefficient  is  independent  of  temperature.  We  have  shown  previously  that 
Si  is  also  a shallow  donor  and  that  if  care  is  not  taken  during  growth.  Si 
will  appear  and  enhance  the  already  considerable  electron  density. 
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In  order  to  intentionally  dope  GaN,  either  to  compensate  the  native 
donor  electrically  or  to  depress  it  via  a chemical  mass  action  law  effect, 
or  to  try  to  grow  p-type  material,  one  would  like  to  identify  the  elements 
that  are  donors  or  acceptors,  shallow  or  deep,  perhaps  their  solubilities, 
diffusion  coefficients,  etc.  But  this  information  cannot  be  obtained  from 
simple  doping  experiments  followed  by  Hall  measurements  because  of  the 
dominance  of  the  native  donor.  From  the  literature,  it  is  known  that  Zn 
and  Mg  can  somehow  compensate  GaN,  but  Zn-doped  material  sometimes  comes 
out  n-type,  sometimes  high  resistivity, with  no  simple  theory  explaining 
the  difference.  Thus  we  wished  to  examine  the  behaviour  of  a number  of 
chemical  elements  simply  to  identify  their  electrical  activity  in  GaN. 

Last  year  we  showed  that  we  could  determine  their  behaviour  even  in  the 
presence  of  the  native  donor  by  using  the  results  of  near  bandgap  low- 
temperature  photoluminescence.  Photo luminescence  is  excited  in  GaN  by  a 
laser  and  the  penetration  depth  of  the  excitation  is  only  v2000A.  Hence 
ion  implantation,  which  deposits  dopants  in  a range  of  this  order  or  some- 
what greater  is  an  ideal  way  to  dope  GaN  for  such  measurements.  Further,  from 
the  fluence  and  the  ion  energy,  the  dopant  profile  is  easily  deduced.  The 
only  problem  remaining  is  to  anneal  out  the  damage  caused  by  the  high  energy 
particles  and  "activate"  the  centers. 

In  pure  as-grown  CVO  GaN,  the  dominant  near  gap  luminescence  is  due  to 
an  exciton  bound  to  the  neutral  shallow  native  donor.  After  implanting  with 
Ar+ this  luminescence  disappeared.  By  studying  the  recovery  of  this  band  vs 
temperature  and  time  of  annealing  (in  flowing  NH^)  it  was  found  that  a 3 
hour  anneal  at  1000°C  reproduced  the  original  intensity.  Furthermore, 
when  other  implants  were  done  producing  other  luminescence  features,  these 
annealing  conditions  invariably  optimized  the  emission  intensity.  Thus, 
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the  annealing  procedure  was  established. 

The  ion  implanted  samples  we  report  here  were  implanted  for  us  by 
Hughes  Research  Laboratories  (with  our  gratitude  to  Dr.  l.  Anderson).  These 
include  C+,  Si+,  Ge+*  Sn+,  Be+,  Mg+,  Zn\  Cd+,  0+,  S+,  Li*,  Na\  h\  N+  and 
Ar  . Several  fluences  were  used  for  each  dopant  spanning  the  10  -10  cm 
predicted  concentration  range  at  the  peak  maximum  and  energies  were  chosen 
to  produce  ranges  of  2000  to  3000  %. 


Low  Temperature,  Near-edge  Photo! uminescence 

Although  our  primary  thrust  this  year  has  been  towards  processing  GaN 
at  high  temperatures  and  high  ^ pressures  in  order  to  gain  some  control 
over  the  ubiquitous  native  shallow  donor,  most  of  the  year  has  been  spent 
setting  up,  testing  and  calibrating  the  new  system  and  data  is  only  now 
beginning  to  be  obtained.  However,  we  were  fortunate  in  that  our  program 
to  characterize  donors  and  acceptors  in  as-grown  CVD  crystals,  in  crystals 
doped  during  growth  and  in  ion  implanted  samples  by  the  use  of  low  tempera- 
ture photo luminescence  spectra  has  provided  us  with  a wealth  of  new  informa- 
tion and  this  in  samples  where  the  native  donor  masks  all  other  impurities, 
at  least  to  electrical  measurements. 

The  spectra  were  recorded  on  a Perkin-Elmer  E-l  f/8  spectrometer  with 

an  RCA  1P28  photomul tipi ier,  using  boxcar  integration  to  improve  the  signal- 

to-noise  ratio.  Excitation  above  the  band  gap  was  by  a pulsed  N,  laser 

, o c 2 

(3371  A,  10  W peak  power,  10  ns  pulses,  100  pps)  focussed  onto  the  sample 

immersed  in  liquid  He.  Estimating  a lifetime  of  ~10”^sec,  we  deduce  a bulk 

excitation  rate  of  4 x 10^  electron-hole  pairs  /cm^/  sec  and  a steady-state 

free  carrier  excitation  density  of  -lO^cm"2. 

J 


4 


The  Interpretation  of  the  results  is  strikingly  similar  to  those  for 
the  II-VI  compound  CdS  and  follows  very  closely  the  pioneering  experiments 
in  CdS  of  Thomas  and  Hopfield,  Wheeler  and  Dimmock,  and  Reynolds  and  co- 
workers.  In  fact  GaN  and  CdS  are  quite  similar.  8oth  are  wurtzite  struc- 
tureswith  high  ionicity  and  both  have  very  similar  electron  and  hole  effec- 
tive masses.  Only  their  bandgaps  are  different. 

Previous  work  has  measured  the  band  gap  of  GaN  by  absorption,  and  two 
of  the  three  expected  free  excitons  by  reflectivity,  and  the  lowest  energy 
free  exciton  has  also  been  observed  in  luminescence.  Infrared  reflectivity 
of  the  coupled  plasma-reststrahlen  system  yielded  an  electron  effective  mass 
of  0.2,  but  a hole  mass  of  about  unity  is  arrived  at  only  by  naive  estimation. 

From  the  bandgap  at  4.2°K  (3.503  eV)  and  the  el ectron  mass  (the  hole 
mass  being  much  larger  is  therefore  unimportant)  the  binding  energy  of  the 
lowest  free  exciton  is  calculated  to  be  30  meV  and  the  measured  value  is  28meV. 
From  the  free  carrier  masses  and  the  static  dielectric  constant  we  expect  the 
ionization  energy  of  the  shallowest  (effective  mass)  donors  to  be  30  meV  and 
the  shallowest  acceptors  to  be  150-300  meV.  We  expect  most  transitions  to 
be  replicated  at  lower  energies  by  the  simultaneous  emission  of  10  phonons 
(90  meV).  We  also  expect  a Haynes  Rule  constant  of  0.1  to  0.2  relating  the 
binding  energy  of  an  exciton  bound  to  an  neutral  donor  or  acceptor  to  the 
ionization  energy  of  that  donor  or  acceptor. 

We  start  our  analysis  by  considering  the  emission  within  300  X of  the 

band  edge  from  a series  of  CVD  GaN  crystals  not  intentionally  doped  and 

grown  under  conditions  where  Si  contamination  is  low,  so  that  the  electrical 
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properties  are  completely  dominated  by  the  native  shallow  donor  (10-10 
cnf3).  More  than  a dozen  emission  bands  are  observed.  Several  of  them  are 
clearly  90  meV  10  phonon  replicas  of  other  bands.  The  highest  energy  band 


at  3.4751  eV  is  present  in  most  samples.  It  coincides  with  the  predicted 
lowest  free  exciton  and  agrees  with  that  interpretation  from  previous  reflec- 
tivity spectra.  Next  come  two  lines  3.4710  and  3.4691  eV  close  to  the 
predicted  position  for  the  decay  of  an  exciton  bound  to  a neutral,  shallow, 
effective  mass  donor.  (These  would  be  called  ^-'nes  1°  CdS).  Two  other, 
but  much  weaker  lines  at  lower  energies, are  associated  with  each  of  these 
two  bound  excitons.  As  in  CdS,  we  assume  that  these  extra  lines  are  due 
to  bound  exciton  decay  at  a neutral  donor,  but  that  the  final  state  of  the  transi- 
tion leaves  the  neutral  donor  in  an  excited  state.  We  note  that  the  extended 
excited  states  of  a donor  are  much  more  effective-mass  like  than  the  compact 
ground  state  which  may  deviate  from  its  effective  mass  description  by  the 
so-called  central -cell  corrections.  Hence,  we  assume  that  the  two  excited 
states  we  observe  for  each  donor  can  be  described  by  a simple  Rydberg  series 
corresponding  to  the  n<*2  and  n»3  states.  From  the  difference  in  energy 

between  these  two  lines  for  each  donor  we  obtain  a value  of  the  electron 

effective  mass  of  0.2,  in  agreement  with  our  predictions.  Finally,  the  dif- 
ference between  the  n*2  and  the  non-hydrogenic  nal  ground  state,  when  added 
to  the  Rydberg  difference  between  n*»andna2  gives  us  the  ionization  energy 
of  each  of  the  donors.  The  results  then  define  the  two  donor  ionization 
energies:  34  meV  (very  close  to  the  expected  ideal  effective  mass  value 

of  30  meV)  and  50  meV.  (We  will  show  below  that  the  first  is  the  native 

shallow  donor  and  the  second  is  Si.)  3oth  of  these  results  produce  a 
consistent  value  of  0.12  for  the  Haynes  Rule  constant  relating  the  exciton 
binding  energy  to  the  donor  ionization  energy. 

Next,  at  lower  energies  we  observe  two  bands  (3.457  and  3.455  eV)  in 
the  vicinity  of  the  predicted  position  for  excitons  bound  to  shallow  neutral 
acceptors.  (These  would  be  1^  lines  in  CdS).  The  intensity  ratio  between 
these  lines  varies  from  sample  to  sample  indicating  that  two  independent 
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acceptors  are  involved,  but  the  lower  energy  line  is  usually  much  stronger. 

No  excited  state  replicas  are  observed.  Hence  we  cannot  deduce  the  ioni- 
zation energies  directly. 

At  yet  lower  energies  we  observe  a very  broad  oscillatory  emission 
complex  which  is  simply  a single  broad  band,  strongly  replicated  by  the  emission 
of  1,  2,  3,  4,  5 and  6 10  phonons.  This  structure  and  its  position  is 
the  signature  of  a donor-acceptor  pair  band.  An  electron  trapped  on  a 

donor  recombines  radiatlvely  with  a hole  on  an  acceptor  yielding  an 

2 

emission  energy  E 3 Eg  - (E^  + E^)  + e /er.  Here  r is  the  separation 

between  the  donor  and  the  acceptor  and  e is  the  static  dielectric  constant. 

From  the  statistics  of  the  pair  distribution  based  on  the  known  native  donor 

concentration  and  from  the  overlap  of  the  donor  and  acceptor  wave  function 

determined  by  the  Bohr  radius  of  the  ground  state  of  the  shallow  native  donor, 

we  estimate  the  peak  of  the  emission  should  correspond  to  a pair  separation 

o 

r between  100  and  200  A. 

Next  we  repeat  the  experiment  at  77°K.  Here  we  expect  the  shallow 
donor  to  be  ionized.  Thus  the  recombination  involves  a free  electron  in 
the  conduction  band  recombining  with  a hole  on  the  acceptor.  Then,  comparing 
the  77°  and  4.2°  results  and  using  the  known  bandgap  variation  with  tempera- 
ture, the  mean  kinetic  energy  of  a free  electron  at  77°  and  an  average  donor 
ionization  energy  between  the  two  shallow  donors  for  E^,  we  deduce  that  the 
acceptor  binding  energy  (the  deeper  acceptor  of  greater  intensity)  is  202  meV 

and  that  the  coulombic  term  in  the  pair  recombination  energy  is  11  t 1 meV, 

o 

corresponding  to  a pair  separation  at  peak  intensity  of  140  A (again  close 
to  the  CdS  value).  From  this  acceptor  energy  and  the  position  of  the  deeper 
acceptor  bound  exclton  line  we  conclude  that  the  Haynes  Rule  constant  for 
acceptors  is  0.1  Hence  the  shallower  acceptor  has  an  ionization  energy  of 
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175  meV.  We  will  show  below  that  the  deeper  acceptor  is  Zn  or  Cd  (both  have 
nearly  identical  binding  energies)  and  the  lower  intensity  acceptor  is  pro- 
bably Ge. 

Note  that  we  now  have  a consistent  picture  establishing  the  Haynes  Rule 
constants  for  both  donors  and  acceptors,  with  the  data  being  internally 
self-consistent  and  within  the  limits  of  the  theoretical  predictions. 

Now  we  turn  to  the  doped  samples.  All  of  these  samples  always  exhibit 
the  3.4710  eV  bound  exclton  corresponding  to  a donor  ionization  energy  of 
34  meV.  This  must  be  the  shallow  native  donor  which  dominates  all  our 
samples  in  Hall  and  resistivity  measurements. 

We  have  previously  shown  that  Si  was  also  a shallow  donor  and  that 
it  was  the  dominant  impurity  in  all  our  early  CVO  grown  crystals.  ThuSj 
comparing  the  spectra  from  these  early  crystals,  with  those  from  the  low 
Si  later  crystals,  and  with  samples  ion  implanted  with  Si,  we  can  identify 
the  second  bound  exciton  peak  at  3.469  eV  with  a Si  donor  whose  ionization 
energy  is  50  meV. 

Continuing  with  the  other  ion  implanted  samples,  one  by  one,  using 
the  Haynes  Rule  constants  already  derived,  and  observing  donor-acceptor 
pair  band  shifts  between  different  acceptors  and  donors  and  high  temperature 
free  electron  transitions,  we  deduce  that  both  Zn  and  Cd  are  shallow  acceptors 
with  identical  binding  energies  of  200  meV.  Ge  is  both  a donor  and  an  acceptor 
with  a donor  energy  of  67  meV  and  an  acceptor  binding  energy  of  175  meV.  S, 

C,  0,  Mg,  8e  are  relatively  deeper  levels  where  Haynes  Rule  should  no  longer 
be  applicable.  We  wish  to  complete  higher  temperature  experiments  (free-to- 
bound  transitions)  before  concluding  whether  these  are  donors  or  acceptors 
and  what  their  binding  energies  are.  Samples  ion  implanted  with  Sn,  Li,  Na, 

H and  N are  also  being  analyzed  currently. 
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Thus,  the  photo  luminescence  data  of  as-grown,  doped  during  growth,  and 
Ion  Implanted  samples,  when  completed,  should  give  us  a comprehensive  picture 
of  the  electrical  behaviour  of  all  simple  doping  elements  in  GaN  in  spite  of 
all  samples  being  electrically  masked  by  the  native  shallow  donor.  In  addi- 
tion, this  method  allows  us  to  easily  monitor  the  native  shallow  donor  and 
may,  via  Zeeman  and  stress  splitting  experiments  shed  more  light  on  the 
physical  structure  of  the  native  defect  itself. 

Electron-Hole  Droplet  Recombination,  Stimulated  Emission  and 
Optically  Pumped  Ultraviolet  GaN  Lasers 


To  investigate  the  excitation  dependence  of  the  photoluminescence 
described  above,  the  nitrogen  laser  was  focussed  to  a very  small  spot  on 
the  samples.  At  these  very  high  pump  densities,  a new  very  broad  emission 
band  was  observed  at  energies  just  below  the  position  of  the  shallow  bound 
exciton  lines.  This  band  was  observed  irrespective  of  the  specific  sample 
doping.  At  the  highest  pump  densities  it  dominated  the  extrinsic  spectra 
and  it  exhibited  line  narrowing  with  increased  excitation,  characteristic 
of  stimulated  emission.  Because  the  emission  consisted  of  both  spontaneous 
and  stimulated  components  occurring  simultaneously  and  making  it  impossible 
to  sort  out  the  spectra,  we  decided  to  study  the  gain  characteristics 
separately.  To  do  this  we  needed  a reasonable  optical  cavity.  As  noted 
earlier,  the  surface  morphology  of  the  standard  CVD  crystals  grown  on  R- 
plane  sapphire  is  not  planar.  Rather  than  remove  the  sapphire  from  such 
samples,  polish,  etch  and  cleave  them,  we  decided  to  use  some  of  the  small 
very  planar  platelets  of  GaN  that  frequently  appear  during  CVD  growth, 
nucleating  on  the  tip  of  the  NH.,  inlet  tube.  In  the  first  experiments  a 
one-sided  cavity  was  formed.  A single  end  mirror  was  fashioned  by  cleaving 
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one  edge  of  the  flat  surface,  or  in  some  cases  a naturally  occuring  facet  served 
just  as  well.  The  exciting  nitrogen  laser  was  focussed  onto  a variable  width 
slit  and  an  image  of  the  slit  was  focussed  onto  the  flat  surface  of  the 
crystal  with  the  length  of  the  slit  perpendicular  to  the  edge  and  the  image 
startling  at  the  edge.  Stimulated  emission,  unencumbered  by  spontaneous 
emission  was  observed  from  the  edge.  8y  varying  the  length  of  the  si  it  the  net 
gain  in  the  cavity  could  be  varied  in  such  a way  as  to  determine  the  gain 
spectrum. 

The  gain  spectrum  and  its  temperature  dependence  are  best  interpreted 

in  terms  of  the  electron-hole  metallic  droplet  recombination  mechanism, 

where  the  very  ionic  nature  of  GaN  enhances  the  stability  of  this  phase 

via  the  interactions  of  the  electrons  and  holes  with  the  very  polar  10 

phonon.  The  large  binding  energy  measured,  10  meV  is  in  reasonable  agreement 

with  a simple  model  for  the  droplet  including  this  interaction  and  results 
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in  a deduced  low  temperature  liquid  density  of  4.4  x 10  pairs  per  cm  and 
1 he  phase  diagram  we  construct  from  the  data  indicates  the  extremely  high 
critical  temperature  of  81 °K,  the  highest  such  temperature  yet  reported  in 
any  semiconductor. 

By  cleaving  a second  facet  and  focussing  the  nitrogen  laser  to  completely 
fill  the  planar  surface  between  the  two  edges,  a normal  optically  pumped 
cavity  is  achieved  and  the  spectra  of  the  light  emitted  from  either  of  the 
two  edges  as  a function  of  pump  density  maps  out  the  Fabry-Perot  cavity 
modes  at  low  pumping  Intensity  and  clearly  indicates  the  onset  of  lasing  in 
one  or  several  of  these  modes  as  the  power  is  increased. 

Laser  threshold  measurements  vs.  cavity  length  indicate  high  over-all 
conversion  efficiencies,  quite  comparable  with  good  GaAs  lasers  at  low 
temperatures  and  losses  in  this  spectral  range  of  only  about  75  cm" ^ . Again, 
when  compared  with  the  normal  losses  of  typical  GaAs  lasers,  about  30  cm’\ 
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it  is  obvious  that  our  GaN  crystals  are  really  very  good,  relatively  defect- 
free  semiconductor  crystals. 

These  results  were  obtained  as  a natural  extension  of  the  photoluminescence 
program  described  in  the  last  section,  but  the  reults  do  not  help  us  very 
much  towards  achieving  the  main  goals  of  this  contract.  Therefore,  we  will 
not  belabor  the  description  of  these  experiments,  the  results  and  their 
interpretation.  However,  we  attach,  as  an  Appendix  to  this  Report,  a 
complete  account  of  this  investigation  for  the  interested  reader. 


APPENDIX 

STIMULATED  BGSSION,  LASER  ACTION  AND  ELECTRON-HOLE  DROPLET  RECOMBINATION 
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STIMULATED  EMISSION  AND  LASER  ACTION 


It  was  reported  that  GaN  has  an  external  photo lumin- 
escence efficiency  of  about  S%  at  78°K  and  O.S%  at  room 
temperature^  . We  have  found  similar  photo  luminescence 
efficiencies  at  4.20-K  (or  77°K)  and  at  room  temperature. 

GaN  has  been  shown,  to  be  a firect  gap  semiconductor  with 
E_  • S.SOO  eV  at  ZaK.  Furthermore,  the  occurence  of  very 
high  gain  (g  » 105  cm'1’)  in  the  stimulated  emission^-  2 ^ 
indicates  that  GaN  is  a potential  candidate  as  a laser 
material . 

r 2 ] 

Stimulated  emission  and  laser  action  have  recently 
been  reported  by  Dingle  et  al.  They  determined  the  excited 
length  and  pump  power  dependence  of  the  stimulated  emission 
by  gain  measurements ^ at  a temperature  of  2°K  only. 

They  did  not  see  any  Fabry- Perot  mode  structure  in  the  laser 
output  and  did  not  determine  the  cavity  lasses. 

The  present  work  represents  the  first  detailed  study 
of*  stimulated  emission  and  laser  action  in  gallium  nitride 
at  both  LHe  (4-.2CK)  and  LN^  (77°K)  temperatures.  In  the 
first  section  the  experimental  arrangements  for  the  optical 
gain  measurements  are  described.  These  have  enabled  us  to 
measure  stimulated  emission  at  various  excitation  lengths  § 
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power  densities.  In  the  last  section  we  describe  the 
results,  including  optical  gain,  the  dependences  of  ex- 
citation length  and  power  density,  threshold  gain,  cavity 
losses,  Fabry-Perot  mode  structures,  and  the  nature  of 
laser  transition  as  due  to  electron-hole  droplet  recom- 
bination. 

A.  Experimental  Arrangements 
We  fabricated.  Fabry-Perot  structures  on  epitaxial 
GaN  and  on  whisker  GaN  crystals  for  the  study  of  stimulated 
emission  and  laser  action.  Since  epitaxial  GaN  is  hard  to 

I 

cleave,  most  good  Fabry-Perot  structures  were  selected 
from  whisker  (needle-like)  GaN  crystals.  The  cavities 
had  lengths  of  400  - L700  urn.  Two  cleaved  R(ll0Z)  plane 
surfaces  defined  the  cavity  length  L,  and  the  excitation  A 
was  on  a (1010)  hexagonal  face  of  the  crystallites.  Such 
a typical  as -grown  structure  is  shown  in  Scanning  Electron 
Microscope  and  optical  microscope  photographs  in  Pig.  l. 
After  organic  cleaning,  the  whisker  GaN  sample  was 
' mounted  on  a sample  holder  with  silicone  grease  and  im- 
mersed in.  the  liquid  helium  dewar.  All  temperatures 

. 

reported  are  sample  temperatures.  For  the  pump  intensity 
used,  sample  heating  is  not  expected  to  be  significant. 

j 
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(b)  Scanning  Electron.  Micrtsccpe  — I 


picture 

Fig.  1 'Die  shape  of  a.  typical  cavity  of  whisker  GaN 
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The  experimental  set  up  for  gain  measurements  of  the 
stimulated  emission  is  shown  in  Fig.  2 . The  pump  source 
was  3371  X light  from  a pulsed  (100  Of  peak  power,  10* 3 sec 
duration  time  at  repetition  rate  of  100  pps)  laser. 

The  cross-section  of  the  laser  beam  is  rectangular  with 

. 

typical  dimensions  of  5 x 0.2  cm  and  the  intensity  dis- 
tribution along  the  long  dimension  is  uniform  to  within 
2%.  The  laser  beam  passed  through  a variable  attenuator- 
neutral  density  filter  (Oriel  Optics  Corp.)  and  was 
focused  to  a line  by  a cylindrical  lens.  An  adjustable 
slit  (Gaertner  Scientific  Corp.)  placed  at  the  focal  plane 
of  the  cylindrical  lens  was  used  to  control  the  length  1 
of  the  exciting,  beam.  The  slit  was  imaged  onto  the  sample 
by  means  of  a spherical  lens . 

Typically  the  excited  region  had  a width  of  less  than 
100  urn  and  a maximum  length  of  about  2.0  mm.  Under  these 
conditions  the  maximum  excitation  intensity  at  the  sample 
surface  was  10  MW/cm  - Considering  the  total  losses  of 
S0%  (include  30%  reflected  by  the  GaN  surface,  10%  collec- 
tive losses  and  30%  transmission  losses  in  the  quart:  dewar 
walls r lenses  and  mirrors) , the  absolute  maximum  excitation 
intensity  onto  the  GaN  sample  surface  was  about  SMW/cm-. 

The  stimulated  emission  emerging  from  the  edge  of  the 
sample,  was  focused  on  the  entrance  slit  of  Perkin-Hlmer 
S-l  monochromator.  The  output  signal  of  the  spectrometer 
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'(a) 


(a) 


Fig.  2 

(a)  Optical  and  (b)  electronic  set  up  for  the  optical 

GAIN  MEASUREMENTS 
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was  detected  by  a.  fast  photomultiplier  C&CA-1P28)  and 
averaged  by  a fast  boxcar  integrator  (PAR  model  160)  with 
an  aperture  time  of  20  nsec,  a time  constant  of  300  nsec, 
and  a time  base  of  2 usee.  The  photomultiplier  should  have 
a risetime  short  compared  to  the  pulse  length  of  10  sec, 
and  should  have  a large  wta-r-i wmw  anode  current  for  ease  of 
subsequent  electronic  signal  processing.  The  RCA-1P23 
photomultiplier  has  the  following  characteristics:  9-stages 
modified  S-S  (104)  response  of  2000-6300  S,  range,  typical 
gain.  o£  2.3  x 10^,  anode  pulse  rise- time  of  2 nsec,  average 
anode  current  of  0.3  mA  and  a dark  current  of  3 nA  at  an 
operating  voltage  of  1000  volts.  Hence,  this  is  a reason- 
ably good  photomultiplier  for  detecting  nsec-pulsed  signals. 
There  are  two  pulse  amplifiers  (Keithley  109)  between  the 
PMT  and  the  3oxcar  to  multiply  the  signal  by  a factor  of 
100.  The  output  signal  was  recorded  by  a L 6 N chart  re- 
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B.  Results 


We  have  observed  stimulated  emission  in  the  0.36  urn 
UV  region  from  GaW  whisker  samples  with  optical  pumping  at 
liquid  He  temperature  as  shown  in  Fig.  3 . Similar  results 
were  obtained  at  liquid  nitrogen  temperature,  as  shown  in 
Fig.  r witlr  the  lines  shifted  about  10  X.  toward  longer 


The  stimulated  emission  spectra  as  a function 


OF  EXCITATION  LENGTH  AT  77°K 


wavelengths..  The  peak,  of  the  stimulated  emission  occurs 
at  longer  wavelengths  than  the  spontaneous  luminescence. 

The  intensity  o£  the  stimulated  emission  increases  more 
rapidly  than  linear  (superlinear)  with  increases  of  either 
the  excitation: Length  or  the  excitation  intensity. 

The  optical  gain  can  be  obtained  from  a plot  of  La  I 
vs.  1.  The  threshold  gain  for  lasing  can  be  determined 
directly  from,  the  plot  of  C vs.  P.  Total  cavity  losses  of 
about  7Q  cm  ^ can  be  estimated  from  the  reflection  losses 
and  internal  losses  as  obtained  from  a plot  of  P^  against 
1/L.  Fabry-Perot  mode  structure  has  been  observed  with 
high  spectral  resolution.  The  nature  of  the  laser  tran- 
sition as  deduced  to  be  by  the  electron-hole  liquid  (EHL) 
model  is  presented  in  the  next  section. 

1.  Excitation  length  dependence  and  optical  gain 

The  stimulated,  emission  spectra  are  shown  in  Figs.  3 
and  4 as  a function  of  excitation  length  at  4.2°X  and  77°X, 
respectively.  Stimulated  emission  can  be  identified  as 
stimulated  by  the-  superlinear  dependence  of  the  output 
intensity  at  long:  wavelengths  with  increasing  excitation 
length.  Only  the  excitation  length  was  varied  to  produce 
the  observed  shift  of  emission  to  longer  wavelengths. 

For  shorter  excitation  lengths,  there  is  a broad 
shoulder  on  the  longer-wavelength  side  of  the  stimulated 
emission.  This  broad  band  emission  which  peaks  at  3390  A 
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and  its  relationship  to  lasing  in  GaN  is  discussed  in  the 
next  chapter. 


Over  the  region  for  which  optical  amplification 
dominates  the  output,  I varies  as  I Cl)  * exp  (gi)  where  g 
is  the  gain  in  cm  ^ and  l is  the  excitation  length.  Using 
this  relationship,  the  gain  can  be  determined  directly 
from  the  slope  of  the  straight  line  portion  (small  l region) 
of  a curve  of  log  I vs.  length  1 as  shown  in  Fig.  5 . it 
indicates  a net  gain  as  high  as  10^  cm”^  at  a wavelength  of 
538S  5 at  the  highest  excitation  power  density  of  about 
S MW/ca^  at  i.2°£.  As  the  excitation  length  increases , 
log  I increases  linearly,  followed  by  a region  in  which  the 
exponential  variation  saturates.. 

The  optical  gain  can  be  seen  to  vary  from  3 x IQ4,  to 
10^  cm’1,  at  1.2  and  from  1 x 10^  to  S x 10*  cm'1,  at 
77°K. 


2.  Excitation  density  dependence  and  threshold  gain 

The  near  bandgap  spontaneous  emission  at  low  excita- 
tion levels  contains  free  exciton  luminescence  and  various 
bound  exciton  lines  in  GaN.  At  high  excitation  levels 
(about  1 MW/ca^)  none  of  these  emission  lines  were  found 
to  be  important  because  they  do  not  lead  directly  to  stim- 
ulated emission.  Instead,  several  new  emission  lines 
dominate  the  near  bandgap  spectra.  They  are  located  lower 
in  energy  than  the  bound  exciton  lines. 
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Detailed,  studies  o£  the  spontaneous  and  stimulated 
emission  spectra  have  been  made  as  a function  of  excita- 
tion level  at  temperatures  of  both  4.2°K  and  77°K,  as 
shown  in  Figs.  6,  7,  and  8.  At  4.2°K,  the  spontaneous 

emission  spectra,  corresponding  to  light  emitted  from  the 
front  surface  of  the  crystal,  exhibits  a sharp  line  at  3372 
& which  has  already  been  identified  as  recombination  of 
bound  excitons.  The  peak  position  of  the  BE  does  not  shift 
with  increasing  pump  power  density.  A braod  shoulder  cen- 
tered about  339S  $ is  presumably  due  to  electron-hole  liquid 
recombination  and  is  probably  the  origin  of  the  laser  tran- 
sition. This  is  discussed  in  the  next  chapter. 

Stimulated  emission  from  the  cleaved  edges  was  mea- 
sured at  right  angles  to  the  pump  beam  and  its  peak  occurs 
at  longer  wavelengths  than  the  spontaneous  emission.  The 
peak  of.  the  stimulated  emission  shifts  to  longer  wavelengths 
with  increasing  excitation  power  density  as  shown  in  Fig. 

9 . These  two  effects  are  characteristic  of  stimulated 
emission  in  a rather  large  class  of  undoped  semiconductors 
such  as  CdS  and  GaAs. 

The  peak  intensities  of  the  spontaneous  emission  at 
3572  & and  the  stimulated  emission  at  3383  $ are  shown  in 
Fig-  10  as  a function  of  excitation  power  density  at  4.2°K. 
The  excitation  dependence  for  spontaneous  emission  ap- 
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PCD  Excitation  Intensity  (relative  units) 

Fig.  10  The  peak  intensities  of  toe  spotanequs  emission  at 
3572$  AND  THE  STTMJLATED  EMISSION  AT  35552  AS  A 
FUNCTION  OF  EXCITATION  POWER  DENSITY  AT  4.2°K 


preaches  a linear  dependence.  The  intensity  of  the 
stimulated,  emission  increases  more  rapidly  than  linear 
(superlinear)  with  increasing  excitation  intensity.  This 
superlinear  dependence  establishes  the  existence  of  stim- 
ulated emission. 

The  peak  half-widths  (FWHM)  of  the  spontaneous  and 
stimulated  emissions  are  shown  in  Fig.  11  as  a function  of 
excitation  power  density  at  4-.2°K.  The  spontaneous  emis- 
sion band  exhibits  a half -width  of  about  40  - SO  X over 
the  entire  range  of  pump-  power  density  described  here.  The 
half-width  of  the  stimulated  emission  peak  decreases 
abruptly  from  about  40  X at  low  excitation  intensity  to 
about  T X at  high  excitation,  intensity  at  the  threshold 
pumping  power  for  laser  oscillation.  Furthermore,  at  low 
pumping  power  there  is  spontaneous  emission  emitted  in  all 
directions,  whereas  at  the  higher  powers  the  emission  is 
stimulated  emission. 

From  the  plot  of  I vs..  P we  can  determine  directly 
the  threshold  pumping  power  density  for  lasing.  The  exci- 
tation dependence  for  different  cavity  lengths,  L * 400- 
1700  urn,  at  both  4.2°K  and  77°K  are  shown  in  Fig.  12  and 


13,  respectively.  The  threshold  pumping  power  densities 
are  t.SS,  1.10,  0.35  and  0.73  MW/ca2  at  4.2°X  and  1.65,  1.10, 
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Q.9G  and  0.80  MV/ car  at  77  K,  corresponding  to  cavity  lengths 
400,  750,  1200,  and  1700  um,  respectively. 
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carrier  absorption,  penetration  losses,  diffraction  and 
defect-center  scattering  losses) , and  the  reflection  losses 
due  to  the  reflectivity  coefficient  R at  the  laser  ends. 

At  threshold,  the  total  gain  in  the  structure  must 
equal  the  sum  of  total  internal  and  end  losses,  i.e., 

threshold,  condition  for  the  laser  cavity 

is  given  by: 


#(g-a-^)  L 9 j. 

-a-L  (ff-a-)L 
e “■  e 1 » 1 


gth  - ax  ^ 1/1  Lnl/R  * 3 Pth 


(10) 


(U) 


Finally,  the  threshold,  pumping  power  density  equation  is 
given  by: 


Pth  « <xj/8  * 1/1  (1/3  Inl/R) 

i 

where  * threshold  gain 

P^  • threshold  for  pumping  power  density 
3 » gain  factor 

• internal  laser  losses 
1 ■ cavity  length 

R • reflectivity  of  the  ends  of  the  cavity 


(12) 
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We  have  studied  the  dependence  of  P ^ on  laser  cavity 
lengthy  L,  as  shown  in  Table  1 . Equation  (125  predicts 
a linear  dependence.  Figure  14  shows  experimental  results 
at  both  4.2°!C  and  77°K.  for  a set  of  GaN  whisker  laser 
cavities  with  different  cavity  lengths  (400  um  to  1700  um) . 

As  can  be  seen,  the  linear  relation  holds  very  well 
at  both  4.2°£  and  77°K.  With  short  pulses  of  10  ns  dura- 
tion and  a repetition  rate  of  80-100  pps,  sample  heating 
could  be  neglected  even  at  the  highest  power  densities. 

This  dependence  indicates  that  the  gain  is  linearly  depen- 
dent on  the  excitation  power  density  both  at  4.2°X  and  at 
77°K. 

The  reflectivity  of  GaN  at  5.43  eV  in  energy  is  about 
IS  % ^ ‘ . Thus  we  can  estimate  reflection  losses  to  be 

about  30  cm"2.  We  note  from  Fig.  14  that,  for  a given 
cavity  length,  e.g.,  4 x 1Q“2  cm,  the  threshold  power 
density  increases  slightly  from  l.SS  MW/ cm2  at  4.2°K  to 
about  1.63  MW/cm2  at  77°K.  The  corresponding  gain  factor 
ft  decreases  from  47.45  cm/MW  to  44.64  cm/MW,  and  the  internal 
loss  factor  increases  front  24.14  cm"2  to  23.39  cm'2.  It 
is  clear  that  the  threshold  pumping  power  densities  for 
laser  action  and  the  total  cavity  losses  are  almost  the 
same  both  at  4.2°K  and  at  77°iC.  P ^ should  increase  at 
temperatures  greater  than  77°K. 
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TABLE  1 


DEPENDENCE  OF 

THE  THRESHOLD 

PUMPING  POWER  DENSITY  P*u  FOR 

LASER  ACTION 

ON  THE  CAVITY 

LENGTH  L 

Ci:  at  4.2°K 

L Cum) 

1/L  CcaT1) 

Pth  Cl) 

PthfMW/ca2) 

1700 

5.9 

IS 

0.7S 

1200 

3.55 

IT 

0.3S 

750 

14.5 

22 

1.10 

400 

25. Q 

51 

l.SS 

CZ]  at  77°  1C 

kJJaL  I/L  Can’1)  pth  ft)  pthQtW/ca2) 


170  Q 

5.9 

16 

0.30 

1200 

3.55 

13 

0.90 

7S0 

14.1 

22 

1. 10 

4 00 

25.0 

35 

1.65 
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Fig-.  14  ExPERtt^fTAL  RESULTS  OF  THE  DEFSIDBCE  OF  THE  THRESHOLD 
POc©  DENSITY  P-,  ON  THE  CAVITY  LENGTH  L AT  BOTH  4.2?K 
AM)  77°K 
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-Perot  mode  structure 


1 


. 


* Calculation  of  mode  separation 

A standing  wave  is  produced  in  a Fabry -Perot 
cavity  only  if  its  length  L contains  an  integral  number 
of  half -wavelengths,  the  wavelength  measured  inside 
the  medium.  Thus,  nL  » m(\/2),  where  X is  the  wave- 
length in  vacuum,  a is  the  index  of  refraction  of  the 
medium,  and  m is  the  integral  order  number.  For  m > 1 
the  wavelength  separation  between  two  modes  is  then 
\^/2L(n  - X(dn/dA)),  where  in  general  n is  a function 
of  wavelength.  For  the  lasing  wavelength  339S  i,  L ■ 

400  um,.  and  n * 2 (assuming  it  to  be  independent  of 
wavelength) , we  can  estimate  a mode  spacing  of  about 
0.8  & 

b . Experimental  results : 

With  high  spectral  resolution  (0.23  $) , the  laser 
lines  do  indeed  exhibit  the  familiar  Fabry -Perot 
cavity  mode  structure,  typical  examples  of  which  are 
shown  in  Figs.  15  and.  16  . We  found  that  it  consists 
of  fourteen  equally  spaced  modes  (mode  *3  and  *13  are 
missing  for  some  unexplained  reason)  whose  separation 
is  about  0.8  • 0.9  in  agreement  with  the  calculated 
value  of  mode  spacing.  The  peak  positions  of  the 
lasing  modes  are  shown  in  Fig.  17  as  a linear  function 
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Fig.  16  Expert-coal  high-resoluticn  spectral  distribution  of 

THE  FaBRY-PeROT  MODE  STRUCTURE  WITH  DIFFERENT  EXCITATION 
POW6?  DENSITY 
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of  the  node  number.  Because  the  cavity  is  not  a 
simple  Fabry-Perot  structure,  there,  are  cavity  length 
L variations  with  different  excitation  penetration 
depths.  This  effect  could  give  rise  to  the  broadened 
mode  structure  we  observed. 


4 


ELECTRON -HOLE  LIQUID  .RECOMBINATION 


A.  Gain  Spectra  And  Nature  Of  The  Laser  Transition 


In  this  chapter,  we  will  try  to  understand  the  nature 

of  the  lasing  transition  in  GaN.  We  have  described  the 

gain  measurements  on  GaN  whisker  samples  in  the  previous 

chapter.  The  origin  of  the  stimulated  emission  is  due  to 

the  broad  emission  band  which  extends  up  to  about  6S  meV 

below  the  free-exciton  line  (E„_)  at  1.2°K  and  the  lower 

gx 

used  excitation  levels.  These  optical  gain  measurements 
show  that  the  emission  spectra  alone  are  not  sufficient 
to  analyte  the  optical  transitions  involved  at  high  exci- 
tation levels.  It  was  found  that  the  optical  gain  in  GaN 
with  its  particular  wavelength  and  intensity  dependence 
can  distort  the  spectra  of  the  emission  peaks  considerably. 
Therefore,  physical  conclusions  as  to  the  transition  pro- 
cesses involved  should  be  based  on  the  spectral  gain  pro- 
file. 

The  gain,  spectra  of  the  stimulated  emission  peak  was 
obtained  by  analyzing  the  luminescence  intensity  vs.  exci- 
tation length  data  as  first  described  by  Shaklee  and 
Leheny'*  J ' . The  complete  gain  spectra  for  GaN  at  both 

and  7T°X  are  presented  in  Fig.  18  as  a function  of 
wavelength  below  the  band  gap. 
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The  gain  distribution  is  asymmetric  with  a long  wave- 
length tail.  Over  the  range  of  pump  power  density  investi- 
gated (0.5  - 5.0  MW/cm2) , there  is  no  measurable  gain  at  or 
above  the  bandgap  for  either  temperature.  In  addition,  the 
peak  of  the  gain  spectra  does  nor  shift  appreciably  when 
the  excitation  intensity  is  varied  by  a factor  of  10.  The 
gain  linewidth,  measured  from  the  high-energy  edge  to  the 
low-energy  edge  of  the  gain  spectra,  is  also  independent  of 
the  pump  power  density  at  1.1°K.  From  the- above  experi- 
mental results*  we  have  strong  evidence  to  presume  that  the 
electron-hole  liquid  (EHL)  exists  in  this  polar  semiconduc- 
tor* GaN-.  This  origin  of  the  laser  transition,  based  on 
recombination  in  a degenerate  electron-hole  liquid  is  a 
more  reasonable  explanation  than  one  based  on  an  exciton- 
exciton  interaction  process  5 ^ . 

R.  Review  of  Electron-Hole  Liquid 
1.  EHL  in  indirect  gan  semiconductors 

The  phenomenon  of  EHL  formation  in  highly  excited 
indirect  semiconductors  like  Go  and  Si  is  now  well  under- 
stood both  experimentally  and  theoretically^  6 ^ . 

It  is  known  that  non- equilibrium  electrons  and  holes 
in  semiconductors  are  bound  in  excitons  at  low  tempera- 
tures by  the  Coulomb  attraction  between  them.  In  fact, 
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almost  all  the  electrons  and  holes  are  bound  into  excitons 

7 

within,  a few  mW/cm  . At  high  exciton  concentrations  ob- 
tained by  increasing  the  intensity  of  the  exciting  light, 
the  interactions  between  them  leads  to  the  formation  of  a 
collective  phase  of  electrons  and  holes. 

In  196&  Keldysh-  first  suggested  that  free  excitons 
(FE)  of  sufficient  densities  could  condense  into  a degen- 
erate electron-hole  plasma,  which  is  electrically  conducting, 
in  order  to  explain  anomalous  photoconductivity  observed  in 
Ge.  Although  the  analogy  is  not  exact,  the  phenomenon  of 
EH1  is  roughly  analogous  to  the  condensation  of  waver  vapor 
into  fog  droplets. 

Exciting  light  incident  on  a pure  crystal  is  absorbed 
near  the  surface,  generating  free  electrons  and  holes; 
these  quickly  form  excitons,  which  diffuse  up  to  a milli- 
meter or  so  away  from  the  surface.  A droplet  embryo  may  be 
formed,  from  a statistical  fluctuation  in  the  exciton  density 
in  a pure  crystal.  If  the  radius  of  the  droplet  embryo  is 
larger  than  the  critical  radius  (-0.QS  um) , it  can  grow  to 
observable  macroscopic  metallic  liquid  drops  containing 
equal  numbers  of  electrons  and  holes. 

Both  the  FE  and  the  EHL  will  decay  due  to  the  anhila- 
tion  of  electron-hole  pairs.  The  photo  luminescence  spec- 
trum of  FE  and  EHL  in  Ge  Cor  Si)  at  an  optical  excitation 
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just  above  threshold  for  the  formation  of  the  drop  phase, 
shows  the  broad  line  shape  characteristic  of  the  metallic 
drops.  Detailed  analysis  yields  the  Fermi  energy  E?  and 
the  condensation  energy  '4. 

The  stability  of  the  EHL  is  strikingly  dependent  on 
band  structure.  In  Ge  and  Si,,  the  complication  in  the 
band  structure  turns  out  to  be  favorable  for  condensation. 
And  this  is  why  the  luminescence  due  to  EHL  was  originally 
observed  only  in  Ge  and  Si.  The  four  equivalent  valleys 
of  Ge  and  six  valleys  of  Si  in  the  conduction  band  effec- 
tively reduce  the  electron  Fermi  energy  by  a factor  of 
(i)*2^  and  (6)~2//2',  respectively.  The  multivalley  struc- 
ture and  mass  anisotropy  of  Ge  and  Si  together  with  their 
indirect  bandgap,  are  responsible  far  the  observed  stability 
of  EHL  in  these  crystals.  Both  Ge  and  Si  are  indirect 
bandgap  semiconductors,  that  is,  electron- hole  radiative 
recombination  must  be  accompanied  by  the  generation  of  a 
lattice  vibration,  a phonon.  This  reduces  the  transition 
probability  and  gives  longer  exciton  lifetimes,  which  is 
favorable  for  drop  formation. 

2.  EHL  in  polar  semiconductors 

For  direct  band  gap  semiconductors  like  GaAs  and  CdS 
the  calculation  of  the  EHL  binding  energies  by  3rinkman 
and  Rice^  was  not  encouraging.  Because  the  small  elec- 
tron effective  mass  and  the  single  conduction  band  minimum 
were  unfavorable  for  EHL  binding. 
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la  the  case  of  GaAs , the  conduction  band  has  a simple, 
spherical,  parabolic  band  shape,  and  the  electron  effective 
mass  is  very  small  («  0.07  mQ) . So  it  can  not  reduce  the 
kinetic  energy  sufficiently,  and  the  condensation  of  the 
non- equilibrium  carriers  in  this  semiconductor  does  not 
lead  to  an  energy  lowering,  at  any  concentration  of  excitons, 
i.e.  the  binding  energy  per  electron-hole  pair  in  the  plasma 
is  close  to  that  of  the  free  excicon.  In  addition,  GaAs 
has  a direct  band  gap  with  no  phonon  assistance  during  re- 
combination. Hence  it  does  not  have  a long  enough  exciton 
lifetime  to  fora  a drop.  So  it  is  hard  to  observe  EHL  in 
GaAs,.  CdS  etc. 

Recently,  it  has  been  shown  theoretically  by  G.  3eni 

r q i 

and  T..  M.  Rice^  J that  the  electron-phonon  interaction 
plays  an  important  role  in  the  stabilization  of  the  EHL 
phase  in  polar  semiconductors.  They  have  calculated  the 
ground- state  energy  of  the  EHL  including  coupling  to  LO 
phonon  , and  found  that  the  polar  coupling  increases  the 
stability  of  the  liquid  phase  relative  to  excitons  . They 
have  also  applied  the  theory  to  Agflr  and  to  direct  bandgap 
semiconductors  like  CdS,  CdSe  and  ZnS;  and  found  that  the 
EHL  is  strongly  bound  in  CdS  and  AgBr  but  only  weakly  so 
in  CdSe. 
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A similar  conclusion,  has  been  reached  by  Keldysh  and 
Silin*-  ^ ^ . Recently,  M.  REsler  and  R.  Zimmermann*-1Q^ 
have  reported  phase  diagram  calculations  of  EHL  in  excited 
direct  gap  semiconductors.  Results  were  given  for  CdS, 

CdSe,  ZnS,  GaAs  and  CdTe.  - 

Experimental  identifications  of  the  EHL  phase  in  such 
highly  excited  polar  semiconductors  as  CdS^11'14-  CdSe*-11^, 

GaAs  15  ^ and  AgBr*-16  ^ have  been  reported  recently  and  these 
measurements  are  in  good  agreement  with,  the  calculations  of 
Beni  § Rice  and  Rosier  §■  Zimmermann  as  shown  in  Table  2 . 

The  dominant  effect  of  both  electron-electron  and 
electron- phonon  interaction  is  to  add  a self- energy  shirt 
to  the  bare  band  energies.  Usually , the  self -energy  is 
screened  by  the  coulomb  interaction.  In  polar  crystals 
screening  is  due  to  lattice  vibrations  as  well  as  to  the 
electronic  excitation.  The  self -energy  shift  which  is  the 

■ 

sum  of  exchange  and  correlation  energies  has  only  a weak 
k-dependence  and  does  not  change  much  with  temperature*-'  ' . 
Therefore,  the  shifted  energy- gap  can  be  calculated  by 
taking  the  lowest -order  diagram  with  the  interaction 
1 screened  by  the  dielectric  function  only  at  zero  tempera- 

ture and  at  the  Fermi -momentum  related  to  the  density  via 
n - kF5/3ir*. 
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I'uruweters  required: 


The  self-energy  shift  due  to  an.  anisotropic  energy 
band  with  longitudinal  and  transverse  effective  masses 
comes  out  to  be  the  polaron  shift  without  excitation. 
The  EHL  energies  should  be  related  to  the  exciton  energy: 

E«  * h * • E*  C13’ 

Finally,  the  binding  energy  of  EHL  with  respect  to  the 
exciton  gas  is  given  by 


_ g . g 

x o,  min.. 


where  EQ  the  minimum  of  the  ground- state  energy  of 

EHL  is  measured,  from  a plot  of  Efl  vs.  n. 

The  standard  way  to  investigate  the  EHL  in  highly- 
excited  semiconductors  is  to  measure  the  luminescence  spec- 
trum. The  EHL  peak  shows  up  as  a broad  band  with  increasing 
excitation.  For  the  direct  band  gap  semiconductors,  gain 
measurements  have  proved  to  be  a very  powerful  tool  for 
investigating  EHL*-17  ^ because  the  EHL  creates  a population 
inversion  resulting  in  large  optical  gain.  The  normal 
recombination  luminescence  spectra  will  be  distorted  due 
to  the  wavelength  dependent  amplification  of  luminescence. 
Therefore,  luminescence  studies  are  unreliable  for  deter- 
mining EHL  properties.  However,  the  gain  spectra  can  be 
used  directly  to  establish  the  distribution  of  excited 
states  in  the  direct  gap  crystals . 
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Since  the  absorption  changes  sign  at  hv  * u(n,T)  for 
all  temperatures,  the  position  of  the  chemical  potential  can 
be  determined.  When  this  cross-over  point  remains  fixed  in 
energy  over  a large  region  of  excitation  intensity  there 
is  strong  evidence  that  EHL  with  an  equilibrium  density 
n(t)  is  formed.  At  rero  temperature,  u(n,T)  for  the  equil- 
ibrium density  n , u , coincides  with  the  minimum  of  the 

oo 

EHL  ground- state  energy  EQ  , and  this  is  the  way  the 

experimental  value  for  the  binding  energy  of  EHL  'l>  is  de- 
duced. 

* - Egx  - uCn0,  0°K)  - A(FE)  - uQ  (15) 

The  line-width  of  the  gain  spectrum  is  usually  taken  as- the 
Fermi  energy  E- , and  the  excited  carrier  density  n is  cal- 
culated from  this  value  for  given  temperature  and  the  known 
band  masses.  Determination  of  the  binding  energy  of  EHL 
and  the  Fermi  energy  from  the  theoretical  prediction  and 
from  the  experimental  gain  spectrum  are  shown  in  Fig.  19  . 

Because  some  parameters  of  GaN  are  unknown  at  present  , 
the  theoretical  predictions  of  the  binding  energy  and  the 
Fermi  energy  for  EHL  formation  in  GaN,'  based  on  the  3eni  5 
Rice  method,  have  not  yet  been  developed. 

The  band  structure  of  GaN  is  very  similar  to  those  of 
CdS  and  ZnS.  So  we  can  estimate  the  values  of  the  binding 


energy  (-7-15  meV) , the  liquid  density  nQ  (-1-5  x IQ13 
cm’3),  the  critical  density  nc  (-0.3 -1.1  x 1013  ca*J),  and 
the  critical  temperature  T.  ('30°K^ 

C.  EHL  Model  Tn  GaN 

Our  data  can  be  best  described  by  an  HHL  bound  by 
10  meV  with  a Fermi  energy  of  55  meV. 

In  Fig. 18  the  energy  difference  between  the  high- 
energy  edge  of  the  gain  spectra  and  the  free-exciton  energy, 
corresponding  to  the  HHL  binding  energy,  is  given  as: 

V (EHL)  » 3.177  eV  (3365  X)  - 3.1 67  eV  (3573  %) 

» 10  meV  3 1.2°K 

At  0°K  the  gain  linewidth  corresponds  to  the  sum  of  the 
carrier  Fermi  energies.  This  is  directly  related  to  the 
carrier  density  nQ  within  the  EHL  drop.  At  high  tempera- 
tures, the  thermal  excitation  of  the  Fermi  system  modifies 
the  distribution  of  carriers,  but  the  gain  linewidth  is 
still  equal  to  the  sum  of  the  carrier  Fermi  energies . The 
gain  spectra  of  GaN  at  77°K  for  different  pump  power  den- 
sities P(%)  are  shown  in  Fig. 20  „ Here  the  gain  linewidths 
are  significantly  reduced  with  decreasing  pump  power 
density,  but  their  half -widths  (FWHM)  of  remain  almost  con- 


stant. 


Gain  (normalized) 


* 


For  GaN , the  combined  Fermi  energy  of  electrons  and 
holes  at  T • 0°K  is 


EyCO)  • Ef  * ^ * i)  (iA,)2'3  (16) 

e h 

* 

where  nr  * 0.2  nr 

0 

■fi  * tV//'“h,x)l/3  * cs  x °-aZ)l/3  - !•<«  ■„ 

Thus, 

SF(0)  - 2. OS  x 10’U  no2/3  (17) 

From  this  equation,  we  can  estimate  the  expected  linewidth 
for  a plasma  having  an  average  density  Nf  corresponding  to 
the  excitation  density,  assuming  the  carrier  lifetime  to 
be  10  10  sec  and  the  excitation  penetration  depth  of  about 
10'S  em- 


it is  clear  that  the  measured  linewidth  at  4.2°X  is 
independent  of  the  pump  intensity,  and  is  much  broader  at 
low  pump  intensity  than  expected  for  a simple  plasma  with 
density  N(cm"3)  as  shown  in  Fig.  21  . The  data  of  linewidth 
SS  meV  are  best  fitted  with  a fixed  density  a * i.l  x 10^3 
cm  3 which  is  the  carrier  density  within  the  drop  at  T * 0°X. 
At  77°K  the  measured  linewidth  drops  rapidly  as  the  pump 
intensity  is  decreased. 

The  Fermi  energy  Ep  is  assumed  to  be  fixed,  consistent 
with  a constant  value  of  nQ-  Strictly  speaking  this  is  true 
only  at  T * 0°K.  For  finite  temperatures , a correction  to 
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the  Fermi  energy  is  necessary  due  to  the  variation  of  E? 
with  temperature  at  constant  density,  and  also  due  to  the 
change  of  n itself  with  temperature  since  n(T)  itself  is 
determined,  by  minimizing  the  total  free  energy  F of  the 
system. 

For  the  case  where  M » nQ  (N  * 10  cm  ) the  pump 
laser  can  generate  an  excitation  density  greater  than  the 
liquid  density  and.  we  observe  less  temperature  variation 
of  the  gain  linewidth.  In  this  case,  the  observed  reduc- 
tion of  linewidth  with  temperature  can  be  determined  by 
calculating  the  thermal  variation  of  the  sum  of  the  elec- 
tron and  hole  Fermi  energies  for  constant  carrier  density, 
i.e.  a(T)  * const.  » nQ.  The  correction  to  Ep  with  temp- 
erature for  constant  density  aQ  is  well-known  from  statis- 
tical mechanics.  This  is  equivalent  to  the  familiar  Som- 
merfeld  term  -ir"(kgT)  /12Ep  for  a free  electron  gas,  and 
is  given  by 


SEp(T) 


v2  cyj* 

11  EfCO)  Ep(Q) 


EpCO) 


(13) 


It  is  shown  by  the  dashed  curve  in  Fig.  22  . 

For  the  case  where  M <a  . the  correction  to  S-  due  to 

O r 

variation  of  n(T)  with  temperature  is  a more  important  cor- 
rection and.  can  be  evaluated  from  thermodynamics,  i.e.,  the 
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Fig.  2z  Variation  of  gain  linbodth- wrm  temperature.  The  circles 

ARE  EXFERIfENTAL  DATA.  THE  DASHED  CURVE  IS  CALCULATES  3Y 
THE  VARIATION  OF  CONSTANT  DENSITY  FERMI  LIQUID,  AND  THE 
SCUD  CURVE  rS  THE  PHASE  DIAGRAM. 


thermal  expansion.  o£  the  e-h  liquid..  We  can  write  an  expan- 
sion to  secondorder  in  T for  the  free  energy  per  pair  F of 
the  EHL  using  the  landau  theory  of  the  Fermi  liquid. 

I 

FCn.T)  *EQ(n)  - %r(a)  (kjT) 2 (19) 

where  £Q(n)  is  the  ground-state  energy  per  pair  as  a func- 
tion of  e-h  pair  density  in  the  drop.  The  second  term 
essentially  corresponds  to  an  entropy  contribution  -TS 
with  the  entropy  linear  in  temperature,  S » yT,  character- 
istic of  a metal.  The  coefficient  y(n)  is  given  by 

TCa)  * Jr  C3nj2/:5(B49  * "W  C20) 

Equation  (19)  can  be  expanded  to  second  order  about  the 
equilibrium  density  nQ, 

F(n,T)  • E0(n0)  * hE"0(n-n0)2  - hy(n)(k3T)2  (21) 

where  primes  denote  derivatives  w.r.t.  n evaluated  at  n0 . 

3y  minimizing  F(a,T)  in  equation  (21): 

ET  * 0 at  ii  * a. 
o o 

5F/5n  - 0 - E”(a-n0)  - %y'(n)(kBT)2  (22) 

The  variation  of  density  with  temperature  is  then: 

a(T)  - n0  * (k^T)  2 (23) 

■ Eo 
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By  combining  Equations  (13)  and.  (23),  it  can  be  seen  that 
the  Fermi  energy  Ep(T)  at  arbitrary  T is 


cm  ^ 9il‘  . 

Ep(T)  - Ef(0)  (1  u «■  5m 


EF  Ei 


(24) 


o o 


For  convenience  in  discussing  the  experimental  data 
we  introduce  the  pre£actors  <5a  and  5^  according  to 


tt(T)  » nQ(l  - BjjdCjT)4) 

EF(T)  - Ef(0)  (1  - dpd^T)2) 


Finally,  we  can  estimate  the  value  of  <3.  and  5r: 

n c 

-2 


(2S) 

(25) 


5a  - 0 . 0145  ± 0.0013  (mev) 


<SH  » 0.0136  T Q.QQQS  (mev) 


-2 


Now,  we  can  estimate  the  critical  temperature  using 
the  equation  as  follows: 


*cCV  *ao  <L  * Wc>  ) 


(27) 


13 


where  nc  is  graphically  estimated  by  the  value  of  1.3  x 10 
cm’3,  which  is  indicated  by  the  arrow-head  in  Fig.  21. 
Therefore,  the  critical  temperature  T.  can  be  determined 
from  Equation  (27)  to  be  31°IC.  The  liquid-gas  coexistence 
curve  in  the  phase  diagram  of  GaN  is  shown  by  the  plot  of 
T/T  vs.  n/nQ  in  Fig.  23. 
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In  conclusion,  by  analyzing  the  spectral  and  temper- 
ature dependence  of  the  laser  transition  we  deduce  that  it 
is  due  to  electron-hole* liquid  recombination.  It  is  iden- 
tified here  for  the  first  time  in  GaN,  a relatively  polar 
semiconductor.  A binding  energy  of  10  meV,  a very  high 
critical  temperature  of  31°K,  and  a low  temperature  liquid 
density  of  4.4.  x IQ3'®  cm”3  are  deduced  and  provide  strong 
evidence  for  the  proposed  mechanism  of  the  transition. 
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